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Abstract 

In two recent articles Salazar and Brenig question the validity of kinetic theory for gran- 
ular gases and fluids, on the based of a supposedly exact hierarchy of coupled equations 
for the velocity moments, which the authors derive from the BBGKY-hierarchy. Their 
derivation contains several errors, which are exposed. Moreover, they support their 
findings with results from direct monte carlo simulation (DSMC) of the Boltzmann 
equation, which supposedly show that Half's homogeneous cooling law only holds for 
times shorter than one mean free time. However, their DSMC results have no physical 
significance as they are carried out at a density, twice close packing. 

The recent comment in |cond-mat/9911276| by Y. Elskens is completely in line with my 
criticism, voiced at the workshop on Granular Gases, that took place in Bad Honnef on 
March 8 - 12, 1999 (see Ref. [|l|). It prompted the publication of the notes, below made at 
that conference, which raise a number of additional objections. 

Salazar and Brenig in Ref. claim to have obtained the following analytic and simula- 
tion results for undriven (freely cooling) granular gases of Inelastic Hard Spheres (IHS): 

(i) A Direct Simulation Monte Carlo (DSMC) solution of the Boltzmann equation, demon- 
strating that the validity of HafF's law is restricted to times less than one mean free 
time to- This law states that the decay of the total energy from an initially homoge- 
neous equilibrium state obeys E{t) = E{0)/{1 + t/t^Y, where = 2dto/{l — a^) in 
a (i-dimensional IHS fluid, and a is the coefficient of normal restitution. 

(ii) An exact analytic relation for the time dependence of certain velocity moments, on the 
basis of which they question the validity of the kinetic theory of granular gases. 
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In the present comment I show that their analytic relation does not hold, and that their 
DSMC results refer to unphysical densities of twice close packing, and consequently have no 
physical relevance. 

Concerning point (i): many authors have shown by means of molecular dynamics ^, 
and of DSMC [^] that Haff 's law is in quantitative agreement with simulation results for 
times up to the clustering transition (density instability), for all gas and fluid densities, and 
for small and moderate inelasticities with restitution coefficient a > 0.7. 

However, the authors of Ref. arrive at quite different conclusions. From their input 
data for solving the two-dimensional Boltzmann equation by DSMC, I extract the following 
values, relevant for our discussion: number of particles, N = 2 x 10^; system size L = lOOcx 
[a is a disk diameter) and mean free path Iq ~ 9a. This yields the unphysical value for the 
coverage or area fraction (p = jirNa"^ / L"^ ^ 1.57, which is about twice (!) the closed packed 
density of hard disks. According to the Enskog theory the value Iq ~ 9cr corresponds to a 
coverage of ~ 0.03, which is not consistent with the authors' data. Moreover Figure 1 and 
the text above Eq. (3) of Ref. p| inform the reader that the mean free time is ~ 350(it, 
where dt is the integration time step in DSMC, and that Haff's law breaks down at time 
300rft ~ (3/7)to. 

We conclude that the orders of magnitude of the simulation results in Ref. have little 
relevance for DSMC and molecular dynamic simulations at realistic values of the system 
parameter. 



Concerning point (ii): their criticism on granular kinetic theory for three-dimensional 
IHS fluid is based on the following incorrect relation for the velocity moments, dtfik{t) = 
— 27rr7,cr^(l — a^)fik+i{t), claimed to be valid in the translationally invariant homogeneous 
cooling state (HCS). The moments are deflned as, 

f^k{t) = (|c/i2|^) = J dxidx2\gu\''ff\xiX2t). (1) 

Here fQ^\xi . . . Xs,t) are s-particle distribution functions which represent the probability 
density for flnding an arbitrary set of s particles in the volume element dxi . . . dxg, where 
xi = {ri,Vi} {i = l,2,...,s) in the HCS. These pdf's only depend on relative distances. 
Furthermore Qij = Vij ■ fij with Vij = Vi — Vj, and fij is a unit vector along the lines of 
centres at arbitrary nonoverlapping positions inside the volume = of the system. The 
moments /ife(^), are intensive quantities, which approach for large V the values {l/Aim'^) x 

/ dvidv2 J dfi2\gi2\'^fo{vi)fo{v2), which are independent of the system size. This holds for 

r (2) 1 

arbitrary long range spatial correlations in /g (12) — /o(l)/o(2) , which decay faster than 
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1/r" with a > 0. The above argument is made to estabhsh that fik{t) is independent of the 
system size. 

To calculate the rate of change dt{A{12)) of an average pair variable (^4(12)), as in 
Eq. (1), I use the second equation of the BBGKY-hierarchy for the reduced distribution 
functions in terms of binary collision operators Tij as derived in Ref. 0, ^ multiply it 
with (y4(12))/A^^, integrate over xi and X2, and express it in the adjoint binary collision 
operator Tij, which reads 

Ti2 = (r^l dcri2\vi2-cri2\6{fi2-CTcri2){b„^^-l). (2) 

Here au = fi2 for two spheres in contact, and h„^^A{riVir2V2) = A{riVir2V2) with 

V* = Vi- |(1 + a){vij ■ crij)aij. (3) 

In fact, the exphcit form of the operator T12 is not needed in the subsequent arguments, but 
only that the operator is non- vanishing when particles 1 and 2 are close together. 
The result for the rate of change dt{A{12)) is then: 

dt{Ai2) = ^ J dxidx2f^\xiX2) {^12 ■ -^r- + Au + 

+ ^ J dXidX2dx3f^^\xiX2X3)Ti3Ai2 

= /C1 + /C12 + /C123. (4) 

Here )Ci and /C12 represent the terms containing respectively ■ d/dfi2 and T12 on the first 
line of (4), and /C123 the term on the second line. Then, taking A{12) = Igul, ICi vanishes, 
as shown in App.A of Ref. [0. The {xi, X2}-integration in /C12 yields a factor of 0{V), as 
5(ri2 — (JCTu) in (2) is constraint to the value |ri2| = cr. Consequently, the contribution /C12, 
introduced in (A7) of Ref. 0. is of 0{1/V), and vanishes for large systems, and so do the 
RHS of Eqs. (A5) and (31) in Ref. H. 

This does not mean that dtfik{t) in Eq. (31) of Ref. should be vanishing. Inspection 
of /C123 in (4) shows first that this term is 0{1) for large V, because the {ri, r2}-integrations 
yield a factor V"^, because T13 poses the constraint {fi^l = a. For an explicit calculation I 
insert (2) in (4), and take the thermodynamic limit. This shows that the contribution of 
/o^^^(123) - /o(2)/o^^^(13)] to /C123 is of 0{l/V). With the help of the collision dynamics (1) 
and (2) I arrive at the explicit result: 

A^i23 = 1 5- / dvidv2dv3 / dfi2dri3fo{v2)fo^\<ycri3,ViV3) 
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I^^i3 ■ ^3! {|^^i2 ■ ri2 - ^(^^13 ■ ^3)5^13 ■ rul'' - \vi2 ■ rul''} ■ (5) 
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The integrations over {viVsfis} are restricted to the precoUision hemisphere, vis ■ fi3 < 0. 
Inspection shows that /C123 is 0{1) and nonvanishing, contrary to what is argued in Appendix 
AofRef. 0. 

In summary, the relation dtHk{t) = —TkfJ'k+iit), for the velocity moments claimed to 
be exact, does not hold as it is based on an incorrect starting point, Eq. (31) of Ref. 0. 
In deriving their Eq.(31) the authors retain a term /C12 which is of 0{1/V). This term 
should be neglected in comparison to dtfik{t), which is 0{1) for large V. Therefore, the main 
simulation and analytic results of this paper are incorrect. 
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